Search for Deconfined Criticality: SU(2) Deja Vu 
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Monte Carlo simulations of the SU(2)-symmetric deconfined critical point action reveal strong 
violations of scale invariance for the deconfinement transition. We find compelling evidence that 
the generic runaway renormalization flow of the gauge coupling is to a weak flrst order transition, 
similar to the case of U(l)xU(l) symmetry. Our results imply that recent numeric studies of the 
Neel antiferromagnet to valence bond solid quantum phase transition in SU(2)-symmetric models 
were not accurate enough in determining the nature of the transition. 
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Within the standard Ginzburg-Landau- Wilson de- 
scription of critical phenomena a direct transition be- 
tween states which break different symmetries is expected 
to be of first-order. The existence of a generic line of de- 
confined critical points (DCP) proposed in Refs. [H, H, [1| 
— an exotic second-order phase transition between two 
competing orders — remains one of the most intriguing 
and controversial topics in the modern theory of phase 
transitions. In particular, the DCP theory makes the 
prediction that certain types of superfluid to solid and 
the Neel antiferromagnet to valence bond solid (VBS) 
quantum phase transitions in 2D lattice systems can be 
continuous. Remarkably, the new criticality is in the 
same universality class as a 3D system of = 2 iden- 
tical complex-valued classical fields coupled to a gauge 
vector field (referred to as the DCP action below). This 
makes the DCP theory relevant also for the superfluid 
to normal liquid transition in symmetric two-component 
superconductors [1]. 

An intrinsic difficulty in understanding properties of 
the A^-componcnt DCP action is its runaway renormal- 
ization flow to strong coupling at large scales and the 
absence of perturbative fixed points for realistic A" 0, Q . 
One may only speculate that the value of A^ might be of 
little importance since the possibility of the continuous 
transition for A^ = 1 is guaranteed by the exact duality 
mapping between the inverted- XY and XY-universality 
classes 0] and for A^ ^ oo it follows from the large- A^ ex- 
pansion for A^ of the order of a hundred. However, there 
are no exact analytic results either showing that in a two- 
component system there exists a generic line of second- 
order phase transitions, or proving that the second-order 
phase transition is fundamentally impossible. The prob- 
lem of deconfined criticaly for the most interesting case of 
A = 2 thus has to be resolved by numerical simulations. 

The initial effort was focused on models of the super- 
fluid to solid quantum phase transitions and U(l)xU(l)- 
symmetric DCP actions [l|, Q . First claims of deconfined 
criticality were confronted with the observation of weak 
first-order transitions in other models Q. While pre- 



senting a particular model featuring a first order phase 
transition does not prove the impossibility of a contin- 
uous DCP yet, it does raise a warning flag. One needs 
to pay special attention to any signatures of violation of 
the scale invariance which may be indicative of a runaway 
flow to a first-order transition even when all other quanti- 
ties appear to change continuously due to limited system 
sizes available in simulations ■ The fiowgram method 
was developed as a generic tool for monitoring such 
runaways fiow to strong coupling and was used to prove 
the generic first-order nature of the deconfinement tran- 
sition in the U(l) xU(l)-symmetric DCP action. A sub- 
sequent refined analysis resulted in the reconsideration of 
the original claims in favor of a discontinuous transition 
for all known models [12, [3 ■ 

Recently the SU(2)-symmetric case has been studied in 
a series of papers fisl . [l6j and an exciting observation 
of a continuous DCP point was reported. However, the 
story seems to repeat itself since renormalization flows for 
the J-Q model studied in Refs. [13, EB] were shown to be 
in violation of scale invariance and, possibly, indicative 
of the first-order transition [l3| . In this Letter we show 
that a runaway flow to strong coupling and a first or- 
der transition is a generic feature of all SU(2)-symmetric 
DCP models analogous to the U(l)xU(l) case 

For our simulations we consider the lattice version of 
the SU(2)-symmetric NCCP^ model [2, d and map it 
onto the two-component J-current model. The DCP 
action for two spinon fields Za, a — 1,2 on a three- 
dimensional simple cubic lattice is defined as 

- - E i(^:,^a,e'^<-> +C.C) 

<ij>,a 



where {ij) runs over nearest neighbor pair of sites z, j, the 
gauge field v4<y> is defined on the bonds, and V x A is 
a short-hand notation for the lattice curl-operator. The 
mapping to the J-current model starts from the parti- 
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tion function Z = J DzDz* DAexp{—S) and a Taylor 
expansion of the exponentials exp{tz*^z^je^^<^^>} and 

exp{tz*jZ^je~*"^<'^> } on all bonds. One can then per- 
form an explicit Gaussian integration over A<y> , Zai and 
arrive at a formulation in terms of integer non-negative 
bond currents J-'^^. We use ^ — ±1, ±2, ±3 to label the 
directions of bonds going out of a given site the corre- 
sponding unit vectors are denoted by /t. These J-currents 
obey the conservation laws: 

4:1 = 0, with 4J ^ jij - 4;),,.^. (2) 

The final expression for the partition function reads 

■Z' = ^ Qsito Qbond exp(-i7j), (3) 
{J} 

where 

H.J = I E ISV..!^ (4) 

a^b: /j,— 1,2,3 

_ TT 

i.a.fi i,^ ' 

The long-range interaction Vij depends on the distance 
rij between the sites i and j. Its Fourier transform is 
given by = 1/ 2 3 ^^'^^(^m/^) ^^'^ implies an 
asymptotic behavior V ^/fij at large distances. 

This formulation allows efficient Monte Carlo simu- 
lations using a worm algorithm for the two-component 
system ^11]. For the flowgram analysis we measure 
the mean square fluctuations of the winding numbers 
(^a.ft) = (Wa.-fi) of the conserved currents /j-'^^, or, 
equivalently, p± = ± ^2,^,?) /L = ((M>|)/L. 

In particular, we focused on the gauge invariant super- 
fluid stiffness, p- measuring the response to a twist of 
the phase of the product zlz2- 

Similar to the U(l)xU(l) case the NCCPi model 
features three phases. Fig. [U characterized by the fol- 
lowing order parameters: 

VBS: an insulator with (zai) = and, accordingly, 
{p+) = = 0. 

2SF: two-component superfluid (2SF) with (zai) ^ 0, 
(p+)^Oand {p-)^0. 

SFS: supesohd (a paired phase [l^) with {zai) = 
0, (2:1^2:2^) 7^ 0, p+ = and p_ ^ 0. 

The point 5 = and t k. 0.468 features a continuous tran- 
sition in the 0(4) universality class. The relevant part of 
the phase diagram is the region of small g close to this 
0(4) point, far away from the bicritical point g^c ~ 2.0 











- 






/ SFS 


- 


2SF 










VBS 






0(4) 




1 , 1 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 

g 



FIG. 1: (Color online) Phase diagram of the SU(2)-symmetric 
DCP action ([l]). I order transitions VBS-2SF are shown as 
solid red line up to the bicritical point (jjc ~ 2.0. 

where SFS phase intervenes between the VBS and 2SF 
phases. The corresponding direct VBS-2SF transition 
has beenproposed to be a deconfined critical line (DCP 
line) @,[|. 

The key idea of the flowgram method [ll[ is to demon- 
strate that the universal large-scale behavior at g — > is 
identical to that at some flnite coupling g = g^aw where 
the nature of the transition can be easily revealed. The 
procedure is as follows: 

(i) Introduce a definition of the critical point for a 
flnite-size system of linear size L consistent with 
the thermodynamic limit and insensitive to the or- 
der of the transition. In our model we used the 
same definition as in Ref. Specifically, for any 
given g and L we adjusted t so that the ratio of sta- 
tistical weights of configurations with and without 
windings was equal to 7.5. 

(ii) At the transition point, calculate a quantity i?(-/j, g') 
that is supposed to be scale-invariant for a contin- 
uous phase transition in question, vanish in one of 
the phases and diverge in the other. Here we con- 
sider R(L,g) = (VF2). 

(iii) Perform a data collapse for flowgrams of R{L,g), 
by rescaling the linear system size, L C{g)L, 
where C{g) is a smooth and monotonically increas- 
ing function of the coupling constant g. In the 
present case we have C{g ^ 0) (x g 

A collapse of the rescaled flows within an interval g € 
[0, 5coii] implies that the type of the transition within the 
interval remains the same, and thus can be inferred by 
dealing with the g = gcoii point only. Since the g 
limit implies large spatial scales, and, therefore, model- 
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FIG. 2: (Color online) Flowgrams for the short-range model. 
The lower horizontal line features the 0(4) universality scaling 
behavior, so that for g < Qc ^ 0.95 all flows are attracted to 
this line. The upper horizontal line is the tricritical separatrix 
(marked as TP). Above it, flows diverge due to the firs-order 
transition detected by the bi-modal distribution of energy. 
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FIG. 4: (Color online) Data collapse for the NCCP^ flows. 
The yellow line is a fit representing the master curve. The 
horizontal axis is the scale reduced variable C{g)L with 
C{g) = (exp(&p) - l)/(exp(&gi) - 1), fe = 2.28 ± 0.02 and 
gi = 1.3. Error bars are shown for all data points. 



independent runaway renormalization flow pattern, the 
conclusions are universal. 

To have a reference comparison, we first simulated a 
short-range analog of the NCCP^ model (|3]) with Vij = 
gSij . The short-range model has a similar phase diagram, 
but with a second order phase transition for small g and a 
first order one at large g. Figure [2] clearly shows that the 
corresponding flowgram cannot be collapsed on a single 
master curve by rescaling the length (shifting the lines 
horizontally in logarithmical scale), and the separatrix 
at the tricritical point (TP) at g w 0.95 is clearly visible. 

Contrary to the short range model we find no such 
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FIG. 3: (Color online) A typical fiowgram of the gauge in- 
variant superfiuid stifi^ness in the NCCP^ model. The inset 
shows a fan of diverging fiows for 0.125 < g < 1.4 



separatrix for the DCP action. As shown in Fig. [3] the 
flows feature a fan of lines diverging with the system size 
and with the slope increasing with g without any sign of 
a TP separatrix. 

One can notice that the NCCP^ flows exhibit a slope 
change, see Fig. |3] (also observed in Ref. for the J- 
Q-model) that might be interpreted as a sign of the evo- 
lution towards a scale invariant behavior (W^) = const, 
possibly achieved at a large enough L. The same feature 
has been observed recently in Ref. [l^, and caused the 
authors to speculate that the NCCP"'^ model features a 
line of continuous transitions for g < 1.25 20]. The cru- 
cial test, then, is to see if the fan of the NCCP^ lines can 
be collapsed on a single master curve (W^) = F{C{g)L), 
where C{g) describes the length-scale renormalization 
set by the coupling constant g. As it turns out, the 
NCCP^ flows collapse perfectly [2l[ in the whole region 
0.125 < g < 1.65 below the bicritical point gbc (see 
Fig. HI). The rescaling function C{g) exhibits a linear be- 
havior C{g) (X g ai small g consistent with the runaway 
flow in the lowest-order renormalization group analysis 
0. This behavior all but rules out the existence of the 
TP on the VBS-2SF hne. 

Though our conclusions directly contradict claims 
made in Refs. [13) 13) 13, the primary data are in agree- 
ment. A data collapse of the flowgram presented in the 
lower panel of Fig. 13 in Ref. Jla| shows the same quali- 
tative behavior as our Fig. [4] [24]. We are also consistent 
with the conclusion reached in Ref. [13] that the slope 
change is an intermediate scale phenomenon and the Neel 
antiferromagnet to VBS transition in the J- Q- model vi- 
olates the scale invariance hypothesis as observed by the 
divergent flow of {Wl). 

The flow collapse within an interval g G [0,5cou] does 
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FIG. 5: (Color online) Evolution towards the bi-modal energy 
distribution with increasing system size indicative of the first- 
order deconfinement transition {g — 1.65). 



not yet imply a first-order transition. What appears to 
be a diverging behavior in Fig. [3] might be just a recon- 
struction of the flow from the 0(4)-universality (at g = 0) 
to a novel DCP-universality at strong coupling. To com- 
plete the proof, we have to determine the nature of the 



transition for g = gcoii- In this parameter range the stan- 
dard technique of detecting discontinuous transitions by 
the bi-modal energy distribution becomes feasible. As 
shown in Fig. O a clear bi-modal distribution develops 
at g = 1.65 which is below the bicritical point gbc and 
within the data collapse interval [0,(7coii]- 

This leaves us with the clear conclusion that the whole 
phase transition line for small g features a generic weak 
first-order transition identical to the one observed in the 
U(l)xU(l) case. Driven by long-range interactions, this 
behavior develops on length scales (x 1/g ^ oo ior small 
g and thus is universal. It cannot be affected by mi- 
croscopic variations of the NCCP^ model suggested in 
Ref. [lal to suppress the paired (molecular) phase. 
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